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^ ■ Abstract 

1^^ , We calculate parameters in the low energy gravitational effective action and the entan- 

^? , glement entropy in a wide class of theories characterized by improved ultraviolet (UV) 

behavior. These include i) local and non-local Lorentz invariant theories in which inverse 
propagator is modified by higher-derivative terms and ii) theories described by non-Lorentz 
^ , invariant Lifshitz type field operators. We demonstrate that the induced cosmological 

j^ . constant, gravitational couplings and the entropy are sensitive to the way the theory is 

modified in UV. For non-Lorentz invariant theories the induced gravitational effective ac- 
tion is of the Horava-Lifshitz type. We show that under certain conditions imposed on the 
dimension of the Lifshitz operator the couplings of the extrinsic curvature terms in the 
effective action are UV finite. Throughout the paper we systematically exploit the heat 
kernel method appropriately generalized for the class of theories under consideration. 
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1 Introduction 

It is believed that the transplanckian physics manifests itself in certain modifications of the 
existing field models that makes the renormalizability of the quantum models more obvious 
if not automatic. String theory which at the present stage is the only consistent quantum 



or such 



theory which unifies gravity with other interactions effectively provides a mechanism 
a modification. Yet, the concrete details of this mechanism remain to be understooq^. One 
possible way of how this mechanism may work is via modifications in the field propagator so 
that, as function of momentum, it decays faster than the standard 1/p^. The loop diagrams 
then will be regularized and one even may hope to get a theory free of UV divergences. In 
this paper we study the consequences of this improved UV behavior of field propagator for 
quantities which appear in the effective gravitational action and determine the gravitational 
interaction as we know it from the every day experience and astronomical observations. These 
quantities are cosmological constant and Newton's constant. More precisely we study how the 
UV divergences of these quantities are modified by the modifications of the propagator. In 
particular, we investigate whether those divergences can be completely removed by appropriate 
choice of the modified propagator. 

For simplicity we consider the case of a scalar field. Generalization to other fields is straight- 
forward. We consider two wide classes of theories. The first class includes local and non-local 
Lorentz invariant theories. The inverse propagator in this case can be represented in the fol- 
lowing general form 

G-\p) = F{p'), (1.1) 

where F{p'^) is some (to be specified) function of momentuno p"^ = p^p^. As a typical example 
of theories of this type we consider a theory with higher derivatives. The inverse propagator 
then is a polynomial of degree k oi p^ , 



2k 



np')=p' + J^- (1.2) 



The theory may include infinite number of derivatives, for example 



„2 



F{p^) =p' + m'e- . (1.3) 

The theories with higher derivatives are usually considered as not well defined due to presence of 
ghosts. In this paper we do not address this issue. We note however that a number of approaches 
have been suggested recently to properly treat higher derivative theories [2], [3], |1], [5]. 



^For gravitational field, however, we note that there exists an alternative point of view recently phrased in [I] 
that pure Einstein Gravity is already UV complete in the way that transplanckian degrees of freedom cannot 
exist. 

^Throughout this paper we work in the Euclidean signature so that p^ > 0. 



The second class of theories are theories with broken Lorentz symmetry. These are Lifshitz 
type theories |6], [7]. The inverse propagator in this case is 

G-\co,p)=co' + Fip'), (1.4) 

where the momentum now sphts on time-hke and space-hke components, p^ = {u, p*) . A typical 
example of function F in this case is 



2k 

Fip') = P' + -4i-u ■ (1-5) 



The second term in (11.51) breaks the Lorentz symmetry for large values of space-like momentum 
p while the symmetry is restored in the regime of small values of p (compared to the mass 
scale m). 

The other quantity which we study in this paper and which is anticipated to be sensitive to 
the UV modifications of the theory is entanglement entropy. Entanglement entropy of a field 
system is defined by tracing over modes that reside inside a surface S. If the total system 
is characterized by a pure quantum state the subsystem inside the surface is described by 
a density matrix p with the von Neumann entropy S = — Tr p In p known as entanglement 
entropy (see |8], [9] recent reviews [10], [H] and references therein). This entropy is non- 
vanishing provided there are short-distance correlations in the system. The presence of such 
correlations has two consequences: a) entanglement entropy is determined by geometry of the 
surface S, to leading order by the area of S and in higher orders by intrinsic and extrinsic 
geometry of S; b) the entropy is UV divergent, in d spacetime dimensions in a theory with 
propagator l/p"^ one has that 

S^N^r^l. (1.6) 

The exact coefficient of proportionality in ( 11. 6p depends on the regularization scheme employed, 
A^ is the number of species. Roughly speaking, the short-distance correlation 1/a'^^'^ between 
2 points separated by interval a in d dimensions is reflected in the 1/e'^"^ divergence of the 
entropy. Entanglement entropy thus provides us with information about the structure of the 
quantum field theory in UV regime [12] . By improving the UV behavior of the theory one would 
get a milder divergence of the entropy. It is believed that in string theory the entropy should 
be completely finite [I3]. One of the goals of this paper is to test this statement on the theories 
with modified UV behavior, (II. ip and (II. 4p . 

Although we are primarily interested in theories with -F(p^) growing faster than p^ our 
results are general and applicable to other cases as well. It should be noted that theories in 
which propagator is a non-trivial function of momentum are not rare. Let us remind the reader 
some of theories of this type: 



i) 4- dimensional brane in spacetime with extra dimensions, for instance in a theory with one 
extra dimension the induced propagator along the brane is described by function p^ 

F(/) = |th(Lp), (1.7) 

where L is the size of extra dimension. 

ii) DGP model [TB], the inverse propagator induced on the brane then can be modeled by 
the function [ 16] 

F{p^)=p'^ + m^ . (1.8) 

iii) non-commutative field theories, 

It should be noted that in the last two examples the theory is modified in infra-red (IR) 
rather than in the UV regime. 

This paper is organized as follows. We first consider the Lorentz invariant theories. In 
Section 2 we compute the vacuum energy in a theory of the type ( 11. ip . This computation is done 
in flat spacetime using the momentum representation of the corresponding heat kernel function. 
Covariant calculations in curved spacetime confirm this result, which should be interpreted as 
induced cosmological constant. Entanglement entropy in Lorentz invariant theories is computed 
in Section 3. This computation is done in fiat space-time as well. The curved spacetime is 
considered in Section 4 where we compute the induced Newton's constant. In Section 5 we 
present general results on the structure of the heat kernel for an operator which is function 
of the Laplace operator. Non-Lorentz invariant theories are considered in Section 6, where we 
first present the general structure of the low-energy gravitational effective action and use the 
scaling arguments to analyze the structure of UV divergences. Then we present a calculation 
of entanglement entropy in this class of theories. We finish Section 6 by giving a precise form 
for the dependence of couplings in the effective action on the UV cutoff. Some technical details 
are collected in Appendices. 

2 Vacuum Energy/Cosmological Constant 

As a warm-up we calculate the vacuum energy density in fiat spacetime in a field theory with 
field operator 

V = F{U), (2.1) 

which is an arbitrary function of the Laplace operator D = —d^d^. In this paper we sys- 
tematically use the method of the heat kernel in order to calculate the effective action. Con- 
sider a quantum bosonic field described by a field operator V so that the partition function is 



Z = det ^' ^ V . Then the effective action defined as 

1 r°° r/s 

W = -- -TtK{s), (2.2) 

2 J,2 s 

where parameter e is an UV cutoff, is expressed in terms of the trace of the heat kernel 
K{s,X,X') =< Xle'^'^lX' >. The latter is defined as a solution to the heat equation 

id, + V)Kis,X,X') = 0, 
K{s=0,X,X') = 6{X,X'). 

In flat spacetime one can use the Fourier transform in order to solve the heat equation. In 
d spacetime dimensions one has 

K{s, X, X') = -i- [ d^p e*PM(^--^'-) e-^^(p') . (2.4) 

In this representation the trace of the heat kernel Tr K{s) = J d'^X K{s, X, X'=X) reduces to 
a single integral 

TrKis) = j^^VP,is), (2.5) 

where V = f d'^X is the volume of the spacetime and we introduced factors 

Pn{s) = -r^ dpp--'e-^^'^'\ (2.6) 

which in the case of Laplace operator V = n and hence F{p^) = p^ give just s~"/^. 
The effective action 

W = -\{e)V, (2.7) 

in this case represents the vacuum energy with the density 



oo 



where e is an UV cutoff which is introduced to regularize the UV divergences. The small s 
behavior of Pn{s) is determined by the behavior of function F{p'^) for large values of p. Clearly, 
the functions Pn{s) are always divergent in the limit of small s no matter how fast function 
FijP) is growing with p. Respectively the vacuum energy A(e) (12.81) is UV divergent for any 
function F{p^). The divergence however can be milder depending on the function F{p^). In 
curved spacetime the quantity A(e) (12. 8 p should be identified with the cosmological constant. 

As an elementary particular case consider the standard field operator F{\3) = D in c? = 4 
for which one gets 

6 



The vacuum energy in this case is determined only by the UV cutoff e. This is standard UV 
divergence of cosmological constant in four spacetime dimensions. 

The UV modification of the field operator however alters the structure of the vacuum energy. 
Consider two demonstrative examples. 

Example 1: The function F{p'^) is polynomial in p^ . 
For simplicity we take d = A and 

F(p2) = p2 ^ //m2 ^ (2.10) 

where m is a new scale that characterizes the UV completion of the theory. In this case we 
have that 



PM = ^ (^1 + v^^e^ ('^^^) - 1) J ' (2.11) 

where $(x) is error function -^ J^ dte~^ |17j . 

In the limit of small s such that m'^s ^ 1 the function P4,{s) behaves as 



^4(^) = ^ f A - \l^ + 2 + 0{mV~s) 1 . (2.12) 



sm? V sm? 
Respectively, the vacuum energy reads 

^'^) = l6(2^(7r-^T-'"'"(T)j' ("^) 

where we omitted terms which are finite in the limit e — )■ 0. The leading divergence is now 
quadratic. There are also linear and logarithmically divergent terms in the vacuum energy. 

As could be anticipated, in the regime sm? ^ 1 and em ^ 1 the function ^4(5) and the 
vacuum energy for the field operator (I2.10p reproduce the result (12. 9 p for the standard field 
operator. 

In the limit of small s in arbitrary dimension d the integral Pd{s) is dominated by large 
values of p. Suppose that for large p the function 

or, in terms of the field operator, V ~ D^/m^^^"^^ . Then the heat kernel e~^^ is invariant 
under the following rescalings 

X^" ^ (5^/^X'' , s^(3'^s. (2.14) 

This symmetry tells us that in d spacetime dimensions 

Trir(s)~— , s^O, 

d(k-l) 

W = -\{e)V, A(e)~^^^;— . (2.15) 



Example 2: F{p'^) is exponentially growing function of p^ . 
We take 

F{p') = m2(eP'/^' - 1) . (2.16) 

It is characterized by two parameters, m and A. Our choice for the function F{p'^) is such 
that the pole of the propagator F~^(p^) is at p^ = 0, so that in the infra-red regime we have a 
massless field equation. 

The integral Pi{s) for the function ( I2.16p is calculated explicitly and is expressed in terms 
of the hypergeometric function. In the limit of small s such that siv? ^ 1 we have 

P4(s) =AM ^ + ^ + 7lnsm2 + J(lnsm2)2 + 0(sln2s) ) , (2.17) 



^12 2 ' 2 

where 7 is the Euler's constant, so that the vacuum energy 

diverges logarithmically in e. For m of order of the cutoff the cosmological constant is de- 
termined by scale A which can be considerably smaller than the cutoff. Thus the problem of 
cosmological constant should be reconsidered in the UV modified theories. 



3 Entanglement Entropy 

3.1 The replica method 

Before proceeding we remind a technical method very useful for calculation of entanglement 
entropy in a field theory described by the Laplace operator. This method is known as the 
replica trick, see ref. [18]. Entanglement entropy is associated with a surface S and is defined 
by tracing the modes of quantum fields residing inside the surface. One considers a quantum field 
ip{X) on (i-dimensional spacetime and chooses the Cartesian coordinates X^ = (r, x,2;*, i = 
l,..,d—2) where r is Euclidean time, such that the surface S is defined by condition x = 
and {z\ i = 1, .., d~2) are the coordinates on S. In the subspace (r, x) it will be convenient to 
choose the polar coordinate system r = rsin(0) and x = rcos(0) where angular coordinate 
changes in the limits < (p < 271 . 

One first defines a vacuum state of the quantum field in question by the path integral over 
a half of the total Euclidean spacetime defined as r < such that the quantum field takes the 



fixed boundary condition ipi^r = 0, x, z) = ■ipo{x, z) on the boundary of tlie lialf-space, 

nM^,z)]= I P^e-^t'^l, (3.1) 

wfiere W[iIj] is tlie action of the field. The surface S in our case is a plane and direction of 
the coordinate x is orthogonal to S . So that the surface S defined by condition x = 0, r = 
naturally separates the hypersurface r = on two parts: x < and x > 0. Respectively the 
boundary data ip{x,z) is separated on ip_{x,z) = ipo^Xjz), x < and 4'+ = i'oi^:^),^ > ^■ 
Tracing over ip-lx, z) defines a reduced density matrix 



p(^i,V^^) = / I)^-V[/(^^,^_)v[/(V;^,^_) , (3.2) 



where the path integral goes over fields defined on the whole Euclidean spacetime except a cut 
(r = 0,x > 0). In the path integral the field ip{X) takes a boundary value ip']^ above the cut 
and il)\ below the cut. The trace of ra-th power of the density matrix (13.21) then is given by the 
Euclidean path integral over fields defined on an n-sheeted covering of the cut spacetime. In the 
polar coordinates (r, 0) the cut corresponds to values (f) = 27rfc, k = 1, 2, .., n. When one passes 
through the cut from one sheet to another the fields are glued analytically. Geometrically this 
n-fold space is a fiat cone C„ with angle deficit 27r(l — n) at the surface S. Thus we have 

Trp" = Z[CJ, (3.3) 

where Z[Cn] is the Euclidean path integral over the n-fold cover of the Euclidean space, i.e. 
over the cone C„. 

Assuming that in (13. 3p one can analytically continue to non-integer values of n one observes 
that — Trplnp = —{ada — 1) In Trp^l^^^i, where p = -^ is the renormalized matrix density. 
Introduce the effective action W[a] = — In Z{a), where Z{a) = Z[Cq\ is the partition function 
of the field system in question on a Euclidean space with conical singularity at the surface S . 
In the polar coordinates (r, cj)) the conical space Ca is defined by making the coordinate 
periodic with the period 2'Ka, where (1 — a) is very small. Thus one has that 

S = {ad^-l)W{a)U=r. (3.4) 

The uniqueness of the analytic continuation of Tr p" to non-integer n may seem not obvious, 
especially if the field system in question is not relativistic so that there is no isometry in the 
polar angle which would allow us without any troubles to glue together pieces of the Euclidean 
space to form a path integral over a conical space Cq, . However, the arguments are given that 
the analytic continuation to non-integer n is in fact unique. These arguments in some form 
have already appeared in the literature [19]. Here we briefly reproduce some version of these 



arguments since this issue becomes important when we consider the quantum fields with broken 
Lorentz symmetry. 

Consider a renormahzed density matrix p = y^ . The eigenvalues of p lie in the interval 
< A < 1 . If this matrix was a finite matrix we could use the triangle inequality to show that 

|Trp"| < |(Trp)"| = 1 if Re{a) > 1 . 

For infinite size matrices the trace is usually infinite so that a regularization is needed to 
regularize it. Suppose e is a such regularization parameter and Tr^ is the regularized trace. 
Then 

|Tr,p"|<l if Re{a)>l. (3.5) 

Thus Tr p" is a bounded function on 

the complex half-plane, Re{a) > 1. In fact, the path integral calculation on a conical space 
(see next section) gives 

Tr,p~e —^^'\ (3.6) 

where S{e) is the entanglement entropy. Clearly 03.61) satisfies (13. 5p . Now suppose we know 
Tr eP°|Q=„ = ZqIu) for integer values of a = n, n = 1, 2, 3, ... Then in the region Re{a) > 1 we 
can represent Z{a) = Tr^p" in the form 

Z{a) = Zo{a) + sm{na)g{a) , (3.7) 

where the function g{a) is analytic (for Re{a) > 1). Since by condition (13.51) the function Z{a) 
is bounded we obtain that in order to compensate the growth of the sinus in (13. 7p for complex 
values of a the function g{a) should satisfy condition 

\g(^a = x + iy)\ <e-'^l^l . (3.8) 

By Carlson's theorem [20] an analytic function bounded in the region Re{a) > 1 and which 
satisfies condition (13. 8p vanishes identically. Thus we conclude that g{a) = and there is only 
one analytic continuation to non-integer n namely the one given by function Zo{a) . 

In order to calculate the effective action iy(a) we use the heat kernel method. In the 
context of the manifolds with conical singularities this method was earlier developed in great 
detail in [21], [22]. In the Lorentz invariant case the heat kernel K{s,(j),(f)') (where we skip 
the coordinates other than angle 0) on regular fiat space depends on the difference (0 — 0'). 
The heat kernel Ka{s,(j),(f)') on space with a conical singularity is then constructed from this 
quantity by applying the Sommerfeld formula [23] 

K^{s, 0, 0') = K{s, (t)-^') + - — / cot —K(s, - 0' + w)dw . (3.9) 

Ana Jy 2a 

10 



The contour F consists of two vertical lines, going from (— vr + ioo) to (— vr — ioo) and from 
(tt — ioo) to (vr + zoo) and intersecting the real axis between the poles of the cot ^ : — 27ra, 
and 0, +27ra respectively. For a = 1 the integrand in (13.91) is a 27r-periodic function and the 
contributions of these two vertical lines cancel each other. Thus, for a small angle deficit the 
contribution of the integral in ( 13. 9p is proportional to (1 — a). 

3.2 The Lorentz invariant theories 

In this section we want to adopt the definition of entanglement entropy to a generic Lorentz 
invariant non-interacting theory described by a field operator P 

V = F(n) , (3.10) 

which can be an arbitrary function of the operator D = —d^dfj,. We have to construct the vac- 
uum state of the theory, an analogue of the state (13.11) . Suppose that operator (I3.10p is a poly- 
nomial of degree N of Laplace operator D . First we note that one needs to specify N functions 
on the boundary at r = 0, ipi{z,x) = ipij = 0, z,x) , ip2{z,x) = dripij = 0, z,x), ..,-?/' at (z,x) = 
d^^Tpij = 0,z,x). This is the right set of data needed to uniquely define the eigenf unctions of 
operator P on a space with boundary at r = . 
Thus, we define the vacuum state by 

^[{^„, a = 1..N}] = I D^{t,z,x) e-s/'/'^^V' . (3.11) 

{ipaiz,x),a=l..N} 

This prescription for the state is consistent with the one introduced by Hawking and Hertog [2]. 
The reduced density matrix is defined as before by tracing over all functions {ipa, a = l..A^} 
defined on the left half, a; < , of the hypersurface r = . The trace of n-th power of the density 
matrix is then given by the unrestricted path integral over fields defined on the n-sheeted cover 
of the Euclidean space, Z[Cn] = det~^''^D. 

In order to calculate the effective action Vr[C„] = —\nZ[Cn] for partition function Z[Cn\ = 
det~^"-F(n) we use the heat kernel. As we have already seen, in fiat Euclidean space the heat 
kernel for operator -F(n) is obtained by the Fourier transform 

K{s,X,X') = ^ /t/^pe^P''(^'*-^"')e-^^(^') . (3.12) 

(27r)'^ J 

Putting 2;* = z'*, i = 1, .., d—2 and choosing in the polar coordinate system (r, 0) that = (p'+w 
we have that p^{X—X')'^ = 2pr sin ^ cos 6 , where p"^ = p^p^ and 6 is the angle between d- 
vectors p'^ and {X'^—X''^). Radial momentum p and angle 6 together with {d—2) other angles 
form a spherical coordinate system in the space of momenta p^ . Thus one has for the integration 
measure / d'^'p = i^d-2 J^ dpp'^^^ J^ dO sin'^^^ 6 where VL4_2 = r((d-i)/2) ^^ ^^^ ^^^^ °^ ^ ^^^^ 

11 



radius sphere in d—1 dimensions. Performing the integration in fl3.12p in this coordinate system 
we find 

^•(-. -.'■) = "g;f(„i!!Ljll,/. fdp/'Ji^ (2rps,u f ) e--(''l . (3.13) 



2 

We used the fact that 



d9 sin'^-^Oe''''''"' = ^/^ {^)^ T{^)J^{(3) . 



For the trace 



we find that 



d'^-'^z I dcj) j drrK{(j) = (j)' + w,s) 
Jo Jo 

TrKis,w) = ^-!^A(E)P,_2(s) , (3.14) 

(47r)2 sm 2" 



where A{T,) = f d'^ "^z is the area of surface S. We used the integra 
when derived f l3.14p . For the integral over contour F in (13. 9p we finqf 



that [21], [22] 



r{^) 



C2 a = ^ / cot — -^^ = -^1 - a^) . 3.15 

Collecting everything together we obtain 

TtK^{s) = —^{aVPdis) + 2naC2ia)A{J:)Pd-2is)) (3.16) 

for the trace of the heat kernel on space with a conical singularity. 

For quadratic differential operator "D = D one has F{p'^) = p^ and P„(s) = s~"/^. Then 
fl3.16p reproduces the well-known result [H], [22] for the Laplace operator 

Tr K^{s) = \ {aV + s 27^aC2{a)A{T) ) , (3.17) 

where V = J drd'^^^x is the volume of spacetime and A(S) = J d'^~'^x is the area of surface S. 
Substituting fl3.16p into equation (12. 2p we obtain that the effective action contains two terms. 
The one proportional to the volume V reproduces the vacuum energy already calculated in the 
previous section. The second term proportional to the area v4(S) is responsible for the entropy. 
Applying the formula (13. 4p we obtain entanglement entropy 

A(S) 



s 



12 ■ {A7c)(d-2)/ 



r"^ ds 

/ -Pd^2is). (3.18) 

Je^ S 



•^ We note that convergence of integration along contour F in complex w plane in the Sommerfeld formula (13.91) 
is not affected by modification of F{p^) compared to standard case F{p^) — p^ . This can be seen by performing 
integration over w of p.lSp prior to integration over radial momentum p. The result of w -integration is then 
universal for all possible F(jP). 

12 



The behavior of Pd-2{s) at small s depends on the behavior of function F{p'^) at large mo- 
mentum p. However, no matter how fast is function F{p^) at large p the function Pd-2{s) is 
always divergent in the limit s — )■ 0. Thus we conclude that no matter how good is the UV 
behavior of the field propagator F'^ij?) the entropy (13 ■4p is always divergent. 

The divergence however can be milder if F{p'^) is growing faster than p^ . We consider two 
particular examples. 

Example 1. FijP') is polynomially growing function at large p, 

np')^^^y p'»^^ (3.19) 



The UV divergence of the entropy in this case is 

*" " 12 ■ (4ir)(''-W2 r(f) e" ' 

This expression has a pole at rf = 2 which indicates that in two dimensions the divergence in 
the entropy is logarithmic. 

The interplay of the standard p^ term and the UV term in the inverse propagator can be 
seen in d = 4 for the function 

F{p')=p' + ^^, (3.21) 



^2A _ ^2 , P_ 

i2 



, , a/TI" '^ '^ ^, i-s. 



The entropy then contains a logarithmic term 



^ = 4^(v^- + m2ln(me) + 0(m2)) (3.22) 

487r ^ e ^ ' 



as well as UV finite terms proportional to m? . 

Example 2. F{p^) is exponentially growing function at large p, 

F{p^) ~ m^eP'/^\ (3.23) 

P2{s) = -Yil + Hsm') + 0{s)), 

where we take d = 4. The entropy in this case has only logarithmic divergences 

S=^^A^\n\em)+'jln{em)). (3.24) 

13 



We conclude that entanglement entropy is always UV divergent no matter how fast the 
function F{p'^) grows. The degree of divergence however changes considerably if one considers 
F{p'^) growing exponentially. This behavior resembles that of cosmological constant studied in 
Section 2. 



4 Newton's constant 

Our calculation so far has been made in flat spacetime. Now we want to make a step forward and 
generalize this to a weakly curved spacetimes. Although it is rather unusual to use a momentum 
representation in curved space we show in Appendix A that it can be efficiently exploited on a 
weakly curved gravitational background so that the first few covariant terms in the heat kernel 
can be reproduced by using this method. 

In Appendix A the trace of the heat kernel of operator F{\I\g) is shown to take the form 



Tr e-^^(°«) = ^^ {P,{s) I d'X^g + P,_,(.) | d'X^g ^i? + ..) 



(4.1) 



where we omit all terms containing more than two derivatives acting on metric. This expression 
generalizes the well-known small s expansion of the heat kernel of the Laplace operator D^ , 

^'"°' = (ii {/"■'xVg + sjdtxVglR^.) . (4.2) 

However, we should note that (14. ip is a decomposition in number of derivatives of metric rather 
than a small s expansion. In the case of quadratic operator D^ the two decompositions, in 
number of derivatives and in powers of s, in fact coincide. 

Thus the effective action W = —\ J^ ^ Tr e~^^^^^^ for a theory described by field operator 
F[ng) to lowest orders in derivative takes the form 

W[g,,] = -A(e) J d'X^ - j^^^^^ / d'X^Rig] , (4.3) 

where the induced cosmological constant A(e) is given by expression (12. 8p and G'iv(e) is the 
induced Newton's constant, 

1 1 1 f^ds ^ , , ,, ,, 

T6;^GM^=(4^12/.T^-^^^^- ^'-'^ 

The Newton's coupling (14. 4p in the effective gravitational action (14. 3 p should be compared to 
the expression (13.180 for entanglement entropy which we have derived earlier in this paper. The 
comparison shows that entanglement entropy ( 13.180 can be expressed in terms of the induced 
Newton's constant (14. 4p as follows 
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This clearly resembles the Bekenstein-Hawking entropy of a black hole even though the entan- 
glement entropy (13.181) can be defined for any surface which should not be necessarily a black 
hole horizon. The fact that there is a precise balance between the UV divergence in entan- 
glement entropy and the UV divergence of Newton's constant in the effective action was first 
proposed in [13] , later this was checked in different approaches [24] , [25] , [26] and even has been 
proven [27] to be a general property of quantum fields. As the present consideration shows 
this balance preserves for a more general theory which is modified in UV by local or non-local 
terms in the action. This in fact was earlier observed in [12] for the entropy of black holes in 
the context of the Kaluza-Klein theories. In these theories the propagator at short distances is 
essentially higher- dimensional that manifests directly in the entropy of microscopic black holes 
which shows a higher- dimensional behavior. 



5 Complete covariant decomposition 

In the previous section we used a non-covariant method and managed to extract several first 
covariant terms in the decomposition of the heat kernel of operator F{ng) . Here we generalize 
this and in fact obtain the complete structure of the heat kernel. 

First, we note that the heat kernel e~°"°s of operator Dg is available in the form of a well- 
known expansion 

^ oo 

Tre--°^ = (^^E-""'''^- (5-1) 

where C„ are the Schwinger-DeWitt coefficients [28j. The first two coefficients are available in 
(14. 2p . In general, the n-th coefficient is an integral of a covariant expression which contains 
n derivatives of metric. Equation (14. 2 p is thus a covariant decomposition of the heat kernel 
in powers of the Riemann tensor and its covariant derivatives. As we show in this section, 
the heat kernel of operator F{ng) can be decomposed in terms of exactly the same covariant 
combinations C„. The coefficients however in this case are not simply the powers s'^~"'/2 b^t 
functions of the proper time s the exact form of which depends on the function F{ng). 
We can formally rewrite the trace of e"*^''^'''' by inserting a delta-function, 



/"OO 

g-sF(n,)^/ dqe-'^^'^Uing-q) (5.2) 

Jo 

and using a well-known Fourier form for the delta-function, we find a very useful representation 

g-.F(n,)^/ dg^-sFi,)^ dae^'^e-°''% (5.3) 

Jo 27rz J_i^ 



"'This representation is valid if there exists one-to-one correspondence between eigenvalues of operators F{Og 



and Dg . For manifolds without boundaries this condition is satisfied. 



15 



which reduces the heat kernel of operator 
F(Dg) to the heat kernel of the well studied 
and much simpler second order differential op- 
erator dg. Contour of integration over a is 
shown in Figure 1. Its form is dictated by 
condition of convergence. Since the expansion 
( 15. ip can be analytically continued to complex 
values of a (provided it has a positive real 
part) we arrive at the following expression for 
the heat kernel 



Figure 1: Contour of a -integration starts at 
negative imaginary infinity, ends at positive 
imaginary infinity and passes zero from the 
positive real side. 



Tre-^^(°«) 



(47r)<^/2 



J2'^n{s)Cn, (5.4) 



n=0 



where 



POO 

r„(s)= / dqe-'^^'^Uniq), (5-5) 

Jo 

M'l) = 7r- / daa^-'/^e'^'^. (5.6) 



Evaluating the a -integral in (15. 6p we should consider two cases: 

i) n > d/2, the corresponding terms in expansion (15. ip are regular terms; 
ii) n < d/2 , the corresponding terms in the expansion ( 15. ip are singular. 

First, we consider the case of regular terms, n > d/2, and additionally assume that d/2 is 
an integer. Then the integral in (15. 6p can be presented in the form 



/ 1 r+ioo \ 

and hence the corresponding integral over q in (15. 5 p reduces tqj 



(5.7) 



(5. 



For n < d/2 the contour in the integral (15. 6p can be deformed to encircle counter clock-wise 
the negative real axis. The r-integration then gives a F-function (see Eq. (8.315) in [T7]). 

d/2-n-l 

The function (15. 5 p in this case is thus represented by the integral 



Tnis) 



l-^ r dqq"'—\-^'^'^^ =P,.,^{S) 
2-n) Jo 



T{d/2-n) 



(5.10) 



'''If d/2 is not integer the g-derivative in (|5.8p should be understood as a fractional derivative. 
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where functions Pk{s) are those defined in fl2.6p . 

Collecting everything together we find that the heat kernel Tr e"'^^*^'^"-' has the following 
covariant decomposition 

^ oo 

Tre--(°^) = ^^-^5:r„(.)a, (5.11) 

where C„ are the Schwinger-DeWitt coefficients for operator D^ [28j and the functions 7^(s) 
are defined as follows 



Pd-2n{s) n < d/2 

'^"^^^ " ^ (_5^)n-'^/2e-«^W n > d/2 

q=0 



(5.12) 



The singular terms n = and n = 1 in equations (15. lip . f l5.12p agree with the expression (14. ip 
obtained in the previous section. 

It should be noted that as can be seen from f l5.12p in general the total number of singular 
terms in the heat kernel and in the corresponding effective action is preserved under a Lorentz 
invariant deformation of operator D^. In fact, an UV divergent term in a theory with operator 
Dg remains to be divergent (although with a different degree of divergence) in a theory with 
operator F{\I\g) . The logarithmically divergent term for n = d/2 in a theory with Og still 
diverges logarithmically in a theory with F{\I\g) . 

Clearly, our result (15. lip . (I5.12p generalizes to the heat kernel of any operator which is a 
function of a differential operator for which the trace of the heat kernel is known in a form 
similar to (15. ip . 

When function -F(p^) is growing not faster than p"^ for large values of p as in (ll.7p . (II. 8p 
we may deform the contours of integration in (15. 3p in a way that it becomes a Laplace type 
transform 

-| /'+200 /"OO 

g-.F(n,)^^/ ^ge-F(«) / rfcre^'^e-^^^ (5.13) 



27ri_ 

The trace of this heat kernel still takes the form (15.111) . (15.121) . Sometimes it may be more 
convenient to use the Laplace type representation since the heat kernel then is expressed in 
terms of heat kernel e^'^'^" where cr is real variable. 



6 Non-Lorentz invariant UV modifications 

As has been understood recently (see [7], [29], [30] and citations to these papers for the latest 
development) the UV modifications should not be necessarily Lorentz invariant. In fact, by 
changing the theory in UV in a non-Lorentz invariant fashion one can get a renormalizable 
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theory which is free of the ghosts. The gravitational degrees of freedom still can be represented 
as a metric 

ds"^ = N^{t, x)dT^ + 2N\t, x)dTdx' + ^^^(r, x)dx'dx^ . (6.1) 

In this paper we use the gauge in which N^{t,x) = 0. The time coordinate r and the spatial 
coordinates {x*, i = 1, .., d—1} then appear in a different way in the field operator 

V = Ar + F(A,,) , (6.2) 

where 

and 

We remind the reader that throughout this paper we work in the Euclidean signature. Operator 
(16. 2 p is a curved spacetime version of the Lifshitz operator [6] used to describe a phase transition 
in condensed matter. 

The quantum theory described by operator (16. 2p is still characterized by an effective action 
which is a functional of the background gravitational field, 

W[N,g,,] = -l r^Tre-^^. (6.3) 

Z J ^2 s 

As in the Lorentz invariant case the effective action can be decomposed according to the number 
of derivatives of the gravitational field. To lowest order, restricting to the number of derivatives 
not higher than two, we gel|j 

W = jdfi (-\{e) + 6i(e)i?('^-i) + h{e)kl + h{t)k^ + h,{e)^^^\ , (6.4) 

where R^'^~^^ is the Ricci scalar for the spatial (rf — l) -metric gij and we introduced the measure 

dii = drd'^-^xN^ (6.5) 

and the extrinsic curvature 

kij = -^9r9ij , k = g'^kij . (6.6) 

The covariant derivative V is determined with respect to metric g^j. The couplings bj, j = 
1, .., 4 are now the gravitational couplings which replace the single coupling, Newton's constant, 
in the Lorentz invariant case. A(e) still has the meaning of the vacuum energy (or cosmological 
constant). In this section we evaluate explicitly A(e) and couplings bj{e), j = 1, .., 4 that depend 
on the UV cutoff e . This calculation is rather involved. That is why it is useful to start with a 
simple consideration which uses some scaling arguments. 



^The other possible terms quadratic in derivatives are integrals of (drN)'^ and KdrN . These terms however 
should be excluded as they are not invariant under the transformation r = t{t'), N{t) — N{t')^ which is a 
symmetry of operator 



6.1 Scaling arguments 

Consider a quantum field described by a field operator of the type fl6.2p where the function F 
is a power of spatial Laplace operator A^, , 

The heat kernel e~*'^ then is invariant under the scaling transformations [31] 

t^Pt, x' ^ /3t^x\ s^(3^s. (6.8) 

Under these transformations the metric invariants transform as 

jl(d-i) _^ p--,Rid-i) , kij -^ (3-^k,^ . (6.9) 

Since the effective action (16.41) is invariant under (16. 8p (provided that e ^ f3e) we immediately 
obtain the dependence on e for the pre-factors that appear in (16. 4p 



XfA - ^ 


bi( 
1 


^^{'^) ' - d+k-l 5 

62(e) ~ 63(e) ~ 


d-l-fe 



d+k-3 y 

e fc 



61(e) ~ 64(e) 

e fc 

(6.10) 

e k 

The dependence on the mass parameter m of the coefficients (I6.10p can be easily restored by 
dimensional analysis. 

We remind that in the Lorentz invariant case there was a precise balance due to the Gauss- 
Codazzi relations, 61 = 63 = —62 . This balance does not survive if the Lorentz symmetry is 
broken as we see from ( ]6.10p . The interesting fact is that the factors a(e), 61(e) and 64(e) are 
always divergent no matter what is the value of the parameter k in (16. 7p . On the other hand, 
the extrinsic curvature terms in the effective action are UV finite if 

k>d-l. (6.11) 

(For k = d — 1 these terms are logarithmically divergent.) In four spacetime dimensions the 
critical value for k is 3. This seems to be consistent with the proposal of Horava [7] for the UV 
modification of gravity action. 

One can analyze the UV divergence of more general terms that contain powers of extrinsic 
curvature and powers of the curvature of (d—l) -metric gij , 



4 /"t//i/C2P7^", (6.12) 



where /C stands for components of extrinsic curvature and 7^ stands for the components of the 
Riemann curvature of metric g^ . The degree of divergence of the term (16.120 is 

d-l + k-2n-2pk , , 

X = — . 6.13 

k 
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If j9 > 1 we find that all terms of the form fl6.12p that contain extrinsic curvature are UV 
finite if k satisfies the condition (16. lip . On the other hand, if fl6.12p contains only the Riemann 
curvature, i.e. p = 0, the term fl6.12p is UV divergent when 

d-l + k 

n < . 

2 

Comparing to the Lorentz invariant case, A; = 1 , we find that the theory with a broken Lorentz 
invariance contains more UV divergent terms. For instance, the term with n = 3 in dimension 
d = 4 is UV finite in the Lorentz invariant case (A; = 1) while it is logarithmically divergent if 
the theory is perturbed by operator (16. 2p with k = 3. 

Thus, we observe an interesting phenomenon: some UV finite terms in the Lorentz invariant 
theory become UV divergent if the theory is perturbed by a non-Lorentz invariant operator. 

Provided condition (16. lip is satisfied we obtain that the gravitational theory in UV contains 
no terms with time derivatives. Up to the second order in derivative the theory in UV regime 
is described by the action 



W, 



UV 



f drd'^-^x^N (-\ + h^R^"-^'^ + &4^^^] • (6.14) 



By the same scaling arguments we find the structure of entanglement entropy in a theory 
with the field operator (16. 2p . The scaling law for the area v4(S) = j d'^~^x is the following 



d-2 . 

A^I3—A. 



Hence we find that ^31j 



S jzr- (6-15) 

e fc 

We immediately note that the predicted UV divergence of the entropy is the same as in the 
Lorentz invariant theory (see Eq. I3.20p . Earlier in [3TJ it was shown that the UV divergences 
of entanglement entropy in a non-relativistic theory described by the Schroedinger type wave 
operator are the same as in the relativistic theory. As the scaling arguments show, this obser- 
vation may be valid also for a more general type of theories with broken Lorentz invariance. In 
the following sections we demonstrate this in an explicit calculation of the entropy. 

As we have seen in the previous sections, in the Lorentz invariant case the entanglement 
entropy is associated with a particular term in the effective action so that the UV divergence of 
the entropy is identical to the UV divergence of the Newton's constant. As a simple comparison 
of (I6.15P and (I6.10p shows the UV divergence (I6.15P is not anymore associated with any term 
in the effective action (16. 4p . Thus the balance between the gravitational action and the entropy 
does not take place if the Lorentz invariance is broken. 
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6.2 Entanglement entropy: non-Lorentz invariant theories 

In this section for simplicity we use a flat space approximation, N{t,x) = 1, gij{T,x) = Sij. So 
that in (16. 2 p two operators A,- = —d^ and A^ = —df commute. As in the Lorentz invariant 
case computing entanglement entropy associated with a surface S we choose {d — 1) spatial 
coordinates {x\ i = 1, ..,d — l} = {x, z"', a = 1, ..,d — 2} , where x is the coordinate orthogonal 
to the surface S and z'^ are the coordinates on the surface S. Then we switch to the polar 
coordinates, r = rsin(0), x = rcos(0). In the Lorentz invariant case the conical space which 
is needed for calculation of entanglement entropy is obtained by making the angular coordinate 
(j) periodic with the period 27ra by applying the Sommerfeld formula (13. 9p to the heat kernel. 

If the Lorentz invariance is broken, as it is for the operator V = At- + F{Ax) , there are 
certain difficulties in applying the method of the conical singularity when one computes the 
entanglement entropy. The difficulties come from the fact that the wave operator V , if written 
in terms of the polar coordinates r and 0, becomes an explicit function of the angular coordinate 
0. As a result of this the operator I) is not invariant under shifts of (p to an arbitrary constant 
(p + w. Only shifts with w = 2Txn where n is integer are allowed. On one hand, this is a good 
news since a conical space with angle deficit 27r(l — n) is exactly what we need to compute 
Tr p" for the reduced density matrix (13. 2p . The bad news however is that in this case one cannot 
apply the Sommerfeld formula since it explicitly uses the symmetry of the differential operator 
under shifts of angle 0. 

The formal solution to the heat equation for operator V with required 27rn periodicity exists 
in the form 

e"'^ -5{r - r) ^ 5(0 - 0' + 2-Knm) . (6.16) 

However, we did not manage to find an explicit coordinate representation for the heat kernel 
(I6.16P which may be useful to our purposes. Even more, we suspect that the action of the 
operator T) on the delta-function in (16.160 may be ill-defined at the location of the tip of the 
cone, r = . 

Since the direct computation is at the moment not available we have to use indirect reasoning 
and scaling arguments. We use a representation similar to (15. 3p . 

/•oo -I rioo 

rj^ -sAr-sF{A^) ^ dq^-sF{g) / da e''" Ti e-^'^^+^^^K (6.17) 

Jo 27rzy_,^ ^ ^ 

The trace of the heat kernel on a conical space with 2?™ periodicity then reduces to the trace of 
exponential of operator (sA,- + ciA^) (we note that a is a complex number). The conjectured 
form for this trace is 

Tr „e-(^^^+'^^-) = \ \^ {a{n)V + b{n)V^A{i:)) , (6.18) 

(47rs)2 (47ra) 2 
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where 

a{n) =n, b{n) = 27^^2(72) , (6.19) 

and C2{n) is defined in (13.15p . 

Below we give several arguments which support our conjecture: 

i) If both s and a are real we can redefine the time coordinate r' = xH-r . Then we have 

' V cr 

that sA,- + (tAj. = (tD, where D is Laplace operator written in terms of coordinates (r',x). 
Denoting by (r', 0') the corresponding polar coordinates we find a relation between two systems 
of polar coordinates 



A — . / "^ QTi /A 'v — r* , / cin /T^ _L nr\ii^ . 



tan </)' = W — tan </) , r =rW — sin + cos^ . (6.20) 

The eigenf unctions of operator D regular at the origin and periodic with period 27rn in angular 
coordinate 0' are easily constructed (see for example [32]), J\m\ (Xr')e~^ . Clearly, these eigen- 
functions if expressed in terms of original polar coordinates (r, 0) remain regular in the origin 
and periodic with period 27rn in angular coordinate 0. The latter follows from the fact that re- 
lation (I6.20p between two angular coordinates is 27r periodic and hence a 2?™ periodic function 
of 0' is automatically a 27rn periodic function of 0. We conclude that the trace of g-C^^T+o-Aa;) 
on a conical space can be identified with the heat kernel trace of operator D on the same space. 
The latter is well-known and given by (I3.17p . Making the coordinate transformation r' = ^J^t 
in fl3Tri) we arrive at flOSj) . fleTTOj) . 

ii) If a is complex the above procedure does not seem to work since r' becomes a complex 
variable so that J\m\ jXr') in the eigenfunction becomes divergent for large values of r for certain 
values of (for instance = vr). However, we note that in the general case of complex s and 
a the operator (sA^ + aA^.) is invariant under independent rescaling of r and x*: 

x' ^ fix' , o- -^ /i^cr . (6.21) 

Computing the trace of g-C^^T+o-Aa,) ^^ space with 27m periodicity we expect that there are two 
contributions to the trace: one from the volume and the other from the area of the surface S. 
Then the general structure of the trace is fixed by the symmetry (I6.2ip and is given by (I6.18P . 
If s = 0" we note that A^ + A^^ is the Laplace operator. The trace of the heat kernel in this 
case is given by (I3.17p . This fixes a{n) and b{n) as in (I6.19p . 



iii) If function F{q) grows slower than q we can use a Laplace transform (I5.13P to express 
g-sF{Aa;) jj-^ terms of e"""^"" , where now a is real variable. The exact form of the Laplace 
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transform can be obtained by deforming the contours of integration over q and a in fl6.17p . In 
this case the arguments of i) apply and we compute the trace of e-C^^^+f^^a:) ^^ ^^ j^_ Deforming 
the contours of integration over q and a back as in fl6.17p we again obtain fl6.18p . f l6.19p . 

We consider this set of indirect arguments as sufficient to support our conjectured form 

dnnD, dnisD. 

Taking (16.181) , f l6.19p and performing integration over a and q in ( 16.171) we find for the trace 
of the heat kernel K{s) = e~^^ (16.171) on a conical space with 27rn periodicity, 

TTKnis) = nTTKn=iis) + j^^27mC2in)A{J:)Pd-2is), (6.22) 

where nTi Kn=i{s) is the bulk contribution. By the arguments presented in Section 3.1 there 
is unique analytic extension of this formula to non- integer n. A simple comparison with the 
surface term in the heat kernel of the Lorentz invariant operator which was obtained in Section 
3 shows that the surface terms of two kernels are identical. We thus conclude that entanglement 
entropy is given by the same formula 

A(E) f^ds 

as in the Lorentz invariant case ( 13.18p . Earlier a similar property of entanglement entropy was 
observed for a non-relativistic theory described by the Schroedinger operator pT]. Here we see 
that this property is valid for a generic theory with the broken Lorentz symmetry. The fact 
that entanglement entropy is the same in a relativistic theory and in a theory with the broken 
Lorentz symmetry, provided the function F is the same in both cases, is surprising. This means 
that the entropy, although sensitive to the UV modifications of the theory, does not distinguish 
whether these modifications are Lorentz invariant or nolo. 



6.3 Gravitational couplings in curved spacetime 

Our method based on equation ( 15. 3p that we have used in fiat spacetime is applicable in the 
case of a "curved space-time", when the lapse function A^ and the {d — 1) -metric Qij are 
arbitrary functions of r and x*, i=l, .. ,d—l. The main technical difficulty, however, is that 
the two operators, A^ and A^, , do not commute in this case and, in particular, g-«('^T+^(^a:)) -/. 
g-sAr .g-sF(Aj;) Qq ^]-^g^^ g^ more general formula, known as the Zassenhaus formula should be used 
(see Appendix C). With this technique applied step by step the calculation of the heat kernel 



^We are aware of the results |33j indicating that entanglement entropy of free fermions may scale faster than 
area law. Apparently this is due to the presence of the Fermi surface in these models. It would be nice to 
understand whether these systems can be studied with our methods. 
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trace Ti K{s) = t^x^-s{At+f(Ax)) jg straightforward but tedious. An essential simplification is 
that we are interested only in the terms which are quadratic in derivatives of the "metric" : the 
lapse function N and metric Qij . Details of the calculation are collected in Appendix C Here 
we just present the final result 



Tr K{s) = I d^ (a^{s) + Si(s)i?(''"^) + B2{s)kl + B^{s)k'' + B^{s 



(VAT)- 
Ar2 



+ .. 



(6.24) 



where djj, = drd'^'^xN ^fg is a covariant measure, R^^~^^ is the curvature for (rf— l)-metric gij, 
V is the covariant derivative defined with respect to g^j , and we skip all terms with more than 
two derivatives. The exact expressions for the functions of proper time are 
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(6.25) 
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- Pd-l{s) - —Pd+3{s] 



-u. (_p, 



d+As, 



(47r)"'/2 
where we used functions Pn(s) defined earlier (12.60 and introduced Pn{s 



PJs) 



r(t) Jo 



dgg5-ie-^^(^) 



Pn(s) 



r(f ) Jo 



dqq^-^F"{q)e 



-sFig) 



(6.26) 



Here we note that in particular case F(Ax) = A^ when Lorentz invariance restores terms Pn(s) 
vanish and terms Pn{s) give needed powers of s that reproduces (14. 2p . 

The UV divergent terms in the effective action (16. 4p are expressed in terms of the functions 
(E25D 



, , , 1 rds ^ , , 

2 J, 2 s 



hk{e) 



ds 



Bk{s) , k = l,2,3,4. 



(6.27) 



As follows from (I6.25P the functions Aq{s) and -Bi(s) are always positive and respectively the 
cosmological constant A(e) and the coupling &i(e) in the effective action (16.41) are negative. If 
the function F{q) is such that F"{q) < then the function B2{s) is negative and respectively 
the coefficient &2(e) in front of kf, in the effective action (16. 4p is positive. The sign of B^{s) 
and respectively h^{t) depends on the form of the function F{q). 
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It is useful to analyze the coefficients 



for function F{q) = q + iv? ^fc^fc 



There are two 

regimes to consider. In the ffist regime s ^ 1/m? . This could be called an infra-red regime 
since in this case one can neglect the higher momentum term in F{q) and approximate it by 
the linear term F{q) ~ g, the resulting theory then is Lorentz invariant and the coefficients 
Aq{s) and Bn{s) are the same as for the relativistic Laplace operator, 

1 , „, 1 .. „ , . „, 1 . 
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(6.28) 



This is the regime where the Lorentz invariance in the effective action (16. 4p restores. 

In the other regime, s ^ 1/m^, one can neglect the linear in q term and approximate 



F{q) 



m 



2-2k^k 



q . This is essentially the UV regime. In this regime we find the approximate 



values for the coefficients Aq{s) , Bn{s 
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(6.29) 



where we introduced Ao{s) = (47r)2s^/^Ao(s) and Bn{s) = (47r)2s^/^i?„(s). 

The two regimes, the infra-red and the UV, are important when we analyze the UV diver- 
gences in the gravitational action (16. 4p . This can be illustrated for the cosmological constant 
A(e) . The UV regularization parameter e, although supposed to be small, can be large or small 
as compared to the mass parameter m which sets the scale at which the UV modification of 
the field operator becomes important. Thus, if e ^ 1/m we are in the regime when the Lorentz 
symmetry is restored and the cosmological constant takes the usual form 

1 1 



Xmie) 



d{A7i) 



(6.30) 



On the other hand, in the UV regime, when e ^ 1/m, the cosmological constant is 

1 k 



Ac/y(e) 
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(6.31) 
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The most interesting behavior is shown by the couphngs &2(e) and &3(e) in front of the 
extrinsic curvature in the effective action (16 .4^ . In the infra-red regime, these couphngs scale as 



^mi^) ~ -dz^ 



1 

and may seem divergent. However, in the UV regime, when e ^ 1/m, the couphng 62 (^3) 
takes a finite value if the UV term in the field operator is characterized by A; > d — 1 . This 
is an interesting example of how some UV divergences can be cured by the higher momentum 
modifications in the propagator. The UV value of the coefficients are calculated explicitly 



uUV _ 
O2 - 


1 ((d^ 
24(47r)'^/2 


-l)k-S 
{k-\ 


^ + 2rf - 3) 


r(^ 


c-d+l\-p( d-2 \ 
2fc-2 J^ y2k-2> 


iJJV _ 

"3 - 


1 {{cP-l)k- 
24(47r)'^/2 {k - 


cP + 'i-d) 


r( 


fc-d+l^-n/ d-2 N 
2fc-2 '^ y2k-2' 


So that the ratio of two couphngs is 
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-d^ + 2d- 

- rf2 + 3 _ 


3 

J' 





(6.32) 



(6.33) 

It is curious that this ratio becomes —1 in the limit of infinite k so that the balance between 
the two extrinsic curvature terms in the effective action (16.41) is the same as in the General 
Relativity. 



7 Conclusions 

In conclusion let us list our main results 

• We have computed the induced cosmological constant and Newton's constant in a Lorentz 
invariant theory with propagator G{p^) = l/FijP'), F{p'^) = p^ + fiv^)-, where /(p^) is 
growing faster than p^ and analyzed the dependence of these constants on the UV cutoff 
e . We have found that generically there are two regimes: in the infra-red regime, when e is 
sufficiently large the induced constants are effectively determined by p^ term in the inverse 
propagator and take the standard form in terms of e. In the UV regime the modification 
/(p2) becomes important and it changes the degree of divergence of the quantities. 

• We have observed that if function /(p^) is exponentially growing, f{p^) ~ m?e^ /^ , the 
cosmological constant scales logarithmically with cutoff, A(e) ~ A.^\v?{me). This behavior 
is radically different from the standard UV behavior of cosmological constant. 



26 



• Although the degree of UV divergence becomes milder when /(p^) is switched on the 
divergence never disappears completely. In fact all terms in the effective action that were 
UV divergent in a theory with propagator 1/p^ remain divergent if propagator is perturbed 
by /(p^). 

• All what we have just said is true for entanglement entropy. The entropy,as we have found, 
is proportional to the area in the same way as in a theory with standard propagator 
although its dependence on the UV cutoff modifies respectively. However, the entropy 
never becomes UV finite even if the modification in the propagator grows exponentially 
with p^ . In fact, as we have observed, the UV divergence of the entropy is the same as 
that of Newton's constant. 

• In the case of theories with broken Lorentz symmetry there is a set of couplings in the 
effective gravitational action corresponding to the geometric terms which are determined 
either by the space-like metric or by the extrinsic curvature. We have obtained the depen- 
dence of the couplings on the UV cutoff. We show that under certain mild restrictions on 
the dimension of the Lifshitz operator which breaks the Lorentz symmetry the couplings 
of the extrinsic curvature in the action become UV finite. This is the first and only ex- 
ample we know when the couplings which were UV divergent in the theory with standard 
propagator become UV finite when the propagator is improved for large momentum. 

• Entanglement entropy in a generic theory with broken Lorentz symmetry described by 
propagator ( 11. 4p has been calculated and we have found that, quite surprisingly, it depends 
on the UV cutoff e in the same way as in Lorentz invariant theory with propagator ( 11. ip 
provided function F is the same in both cases. The area term in the entropy thus does not 
seem to be sensitive to the nature, Lorentz invariant or not, of the UV modification of the 
theory. Also we note that propagators of the type (11. 4p are typical for condensed matter 
systems. We anticipate that our results of Section 6.2 may have important applications 
in studying entanglement entropy in these systems. 

Our findings may have interesting phenomenological applications which we plan to discuss 
elsewhere. 
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A Heat kernel of operator F{\I\g) in momentum repre- 
sentation 

In this appendix we perform perturbative calculation of the curvature term in heat kernel. First 
we decompose the (i-metric 

where S"'^ is metric on flat spacetime and h°'^ is a small perturbation. Then the covariant 
operator 

is presented as a sum of the Laplace operator in flat space and a perturbation 

6n = -d^{h''^{X)dp) + i(9"/i(X))9„ , (A.l) 

where h{X) = 6ai3h°'^{X) . In what follows it is useful to apply a Fourier decomposition for the 
perturbation 

r/3(X) = — ^ f d'^k e-'^^h^^ik) , (A.2) 

where kX = k^Xf^ . We will assume that the function F{n) can be presented as a Taylor series 

F(n,) = Y: ^^O- . (A.3) 

n=0 

A variation of n-th power of the Laplace operator is 

ra-l 
j=0 

Respectively we have for the variation of F{\I\g) , 

n=0 ■ j=0 

To flrst order in perturbation we can represent the heat kernel in the form 

K{s, X, X') = Ko{s, X, X') + K^{s, X, X') , (A.6) 

where 



Ko{s,X,X') = j 



e-^vix-x') ^-sFip-) (A^7) 



d P ^-ip(X-X') ^-sF{p^) 

(27r) 
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is the heat kernel in flat space, {dg + \I\q)Kq{s,X,X') = 0, and the term 

K^{s,X,X') = I ^^e-^P(^-''')K,{s,p,X'), (A.8) 

is due to the perturbation. The heat equation to first order in perturbation reduces to equation 

{ds + F{no))Ki + 5F{U,)Ko = (A.9) 

which can be easily solved to give the function Ki{s,X,X') in the form ( 1A.8I) with 

K,{s,p,X') = I 0-^e-'''^'D{p,kmp,k)H{p,k) , (A.IO) 

D{p,k) 



Q-sF(p^) _ ^-sF{(p-kf) 



F{p^) - F{{p - kY) ' 

n=0 ' j=0 

H{p, k) = {p^{p - k)p - ^<5„^r (p - k)„)K'^{k) . 
The sum over j in expression for function T{p, k) can be easily evaluated 

-fi^'i(n-j~l)((^_ ^^2y _ ((p-fc) )"- {P T 
j=0 

Noting the Taylor series 



B/)''"-j-'((p-fc)y= "p-*y'::j^^ 






n=0 






n=0 



we obtain a simple expression 



^(,..). ^«p^);)-y) (A.n, 



so that the product of functions D{p, k) and J^{p, k) in (lA.lOp becomes 

-sF{p^) _ -sF{{p-kf) 

Dip, k)HP, k) = ^, _ ^^ _ ^^, . (A.12) 

The expression (1A.10I) with ( 1A.12I) is linear in perturbation h'^^ . On the other hand, it is 
nonperturbative in momentum A;" or, in other words, in derivative of /i°^(X). Using (lA.lOp 



we can in principle obtain the complete result for the part of the effective action that is linear 
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in the perturbation. Written in a covariant form this would give us all terms in the effective 
action which are linear in the Riemann tensor and containing arbitrary number of derivatives. 
Since we present the effective action in a form of decomposition in number of derivatives of 
metric the expression for the heat kernel Ki{s, X, X') flA.10|) should be decomposed in a Taylor 
series in momentum fc" . First we note that 

sF(p^) _ sF{(p-kf) 

+ la\[e-''''\{p-kf-p^) + i<?,|e-"'"''|((p- *)= -p=)' + .. . 
With this formula we find a decomposition in k" 

{D-J^- H){p, k) = A(°)(p) + A«(P)^° + ^S(p)^"fc'' + •• , (A.14) 

where 

A^'\p) = d,2[e-'^^p'^]p^pph^^{k), (A.15) 

AW(p) = dAe-'^^'^-p^hiik) - h{k)p^) - aj.[e-^^(^')]p„p^p,/.^^(fc), (A.16) 

41(P) =ldl.[e-'^'~^\6^pp,p^h'"^{k)^ Kk)p^pp + 2p^p„h"p{k)) (A.17) 

-r\dl,[e-^^^^\^PpP,p^W'^{k). 

Now we are in a position to calculate the trace 

TtK,{s) = ld''XK,{s,X,X'=X) = jd'X I -0^^K,{s,p,X). 

Using the formulae for the integration over p which we list in Appendix B and returning to the 
coordinate representation we find 

In terms of perturbation over flat metric the Ricci scalar reads 

R = -d^dph"^{X)+d''d^h{X) . 

Thus the last term in (lA.lSp is proportional to the Ricci scalar. Taking that in this approxima- 
tion in /i"^ we have that ^/g = 1 — lh{X) + 0{h'^) and R = y/gR + 0{h^) we finally obtain 



Tri^(s)= TTKo{s)+TTK,{s) = -^^(^Pd{s) jd''X^ + Pa-2{s) jd'^X^^R 
for the trace of the heat kernel of operator FiOg) . 
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B Some Integrals in Momentum Space 

For d > 2 the following is true 



j^,ld'pdAe-'^^\.P, = -J^J-fP^is), (B.2) 

jly j d'pdl^e-^'^^'^p^p, = J^J-fP^-^is), (B.3) 

j^j d'pd'^.le ^^(^ ^]PaP,P.P, = - ^^^y/2 Y —P'i-^i-')^ (B-4) 

j^, J d'pd;,[e-^^^^\^,p^,..p^,^^, =0, n,meZ, (B.5) 

where O Pn(s) = y^ /q" c^pp^-^ e-^^(p'). 

C Heat kernel of operator A^ + F{Ax) 

In this appendix we want to comment on explicit calculation of the effective action on a "curved 
spacetime", in this case the lapse function N and the (rf—l) -metric gij in (16. ip are arbitrary 
functions of r and a;*, i = 1,2, .., d—1. For the heat kernel of operator At-+F{Ax) we want to 
use a representation similar to (16 .17^ . Then by introducing the constantly rescaled lapse function 
N'^{t,x) = -N'^{t,x) the heat kernel e-(«^T+o-Aa:) ^^^j-^g jj-^^^q ^^j^g heat kernel e"""^ of covariant 
(i-dimensional Laplace operatoio □ = A^ + A^. built from a new metric G^^ = {N'^,gij}. The 
advantage of this representation is that we can use the well-known Schwinger-DeWitt expansion 
of the heat kernel of the Laplace operator. 

The formula (I6.17P however cannot be directly applied since the operators A^ and A^, do 
not commute if the lapse function A^ and the metric gij are nontrivial functions. In this case 
one has to decouple exponents of these operators, i.e. represent Q-'^i^^+^i^^)) as a product 
of two exponents e"'*'^^^^-' and e"**^^ , then use the representation (I6.17P for e"*'^*^^"'-' and 
finally combine two exponents e~*^^ and e"""^"" into one ^-i^^T+trA^:) _ r^^^^ noncommutativity 
of A^ and A^, at each step leads to pre-factors which are local differential operators and are 
constructed from commutators of the operators. Here one can use a Zassenhaus formula. The 
Zassenhaus formula is analogue of Baker- Campbell-Hausdorff formula which arranges the result 
in a different way suitable for our needs 



^ We note that in contrast to the case of a cone space for smooth manifold case the heat kernel of Lapla- 
cian with complex metric is defined straightforwardly. And all further calculations are also justified since the 
underlying DeWitt's procedure of obtaining heat kernel expansion does not rely on realness of metric coefficients. 
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where A,B- some operators and C„ are linear combinations of n commutators of operators 
A and B, for instance one has [3l] 



C2 = -l[A,B] , Cs = k2[B, [A,B]] + [A, [A,B]]) . 
2 o 

In our case we use both direct and transposed variants of the Zassenhaus formula and 
expand exponents with commutators (with C2, C3,.. ) in series. This allows us to effectively 
decouple and recouple exponents of noncommuting operators in question. Here we choose to 
place prefactors to the left of heat kernels assuming all derivatives act to the right 

/»00 /"ZOO 7 



where 



Q{s, a,q) =1 - [A., A.] l^-s'F'{q) - -sa) + [A., [A., A,]] [-s'F'{q) - -s'a] 
+ [A,., A.][A., A.](^.^F'(g)F'(g) + ^s'a' - ^sVF'(g) - ls'F"iq) 
+ [A., [A., A.]] ( - ls'F"iq) + ^s'F'iq)F'iq) - ^sa') +Oi[,f). (C.2) 

We note that our aim is to present trace of the heat kernel (IC.ip as an expansion in number of 
derivatives acting on the metric components, i.e. on functions N{t,x) and gij{T,x). Commu- 
tators in (1C.2I) are differential operators whose coefficients are (products of) derivatives of the 
metric. Counting number of derivatives in these coefficients one observes that the commutator 
[A^., A^] is proportional to at least one derivative which acts on metric components. Each ad- 
ditional commutation with either A^^^ or A,- increases the number of such derivatives. In this 
context 0([ , ]^) produce terms which contain at least three derivatives. Such counting is closely 
related to the notion of background dimensionality [5^I5B]. 

When being integrated f lC.l|) not all terms in (1C.2P are independent which makes such 



representation of Q{s, a, q) nonunique. To reduce to unique representation we use integration 
by parts over q to remove all explicit dependence on a which enters (1C.2I) in positive powerqj. 
The reduced form of operator Q then reads 

Q{s,q) = 1 - ^[A,,A,][A,,A,]s'F'\q) + 1[A., [A,,A,]]s'F'\q) + 0{[,f). (C.3) 



^Positive powers of a in (jC.21) can be represented as powers of dq acting on e''^ which after integration by 
parts effectively act on e^^^'^''^ leading to powers of s and derivatives of F[q) . Note also that no terms at g = 
arise while integrating by parts that can be seen from the precise representation of (jC.ip in terms of contour 
integrals. In context of the formal derivation given in Section 5, these boundary terms (after taking the trace) 
will be proportional to integrals over a of sufficient negative powers of the latter. These integrals vanish since 
they start and end at infinities and does not pass over zero or a cut in complex cr -plain. 
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In particular, in the reduced representation (IC.3P all terms containing F'{q) disappear. This is 
just what one expects since in the Lorentz invariant case (when F'{q) = 1, F"{q) = F"'{q) = 
.. = 0) the operator Q = 1 as can be seen from flC.l|) . 



Commutators of operators A^ and A^; are calculated straightforwardly: 

[A,, A.] [A,, A,] = l6^fc*^fc""V.V,V;V„92 - 32^F^V'lniV V^V.-Vz^^ (C.4) 



1 



+16^^ Vnn A^V-'ln A^ ViVjd^ + 0{Tf 



[A,, [A,, A.]] = -8A;*^ fc'"* V, V, Vz V^ - 8^ V"F^' V, V, V^9, + 24^^^ V"ln N V, V, V„a, 

+8^V'"VnnAr ViV^^^ - le-^V^ln ATVnn AT ViV^92 + 0{DJI^). 

]\^ TV"' 

In expressions above kij is the extrinsic curvature of constant time hypersurfaces fl6.6p . V is 
space-like covariant derivative and 0{DJI'^) denotes terms - differential operators - which in its 
coefficients contain at least 3 derivatives acting on the d-dimensional metric. These terms do 
not affect the result for heat kernel trace which we restrict here to 0{DJI^) order. 

As we have already noticed the operator Q-i'^'^r+a-A^) ^^^^^ -j^^ viewed as the heat kernel e~°'^ 
of Laplace operator D for the new metric G^,, = {N'^,gij} with the rescaled lapse function 
N'^{t,x) = -N'^{t,x). The small a expansion of this heat kernel is well-known [28] and can be 
used in our calculation, 

e d{T,x]T,x) = — -—j^ e 2. 2_^a a.n,{T,x;T ,x ) , (C.5) 



[ATTay/^ 



ra=0 



where a(r, x; r',x') and D{t,x;t',x') are world function and van Vleck-Morette determinant 
(computed for metric Gf^i, ) respectively [28] . The first two Schwinger-DeWitt coefficients dn at 
coinciding arguments read 

ao(r, x;t,x) = 1 , Si (r, x; r, x) = -R{G) = I (R'^'\g) + P - kfk'' + t.d.) , (C.6) 

00 

where t.d. stands for total derivative terms which we neglect hereafter. 

The next step is to calculate the trace of Q (lC.3p . (lC.4p acting on expansion (IC.Sp . Restric- 
tion to terms which contain not more than two derivatives of the metric considerably simplifies 
the calculation. 

The first term in Q (IC.3P acting on fIC.SP gives 



^^^|d/i[l + al{R^'-'\g) + ~k^-~kik-) 

where dfl = drd'^^h: N ^/g . 
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(C.7) 



The rest two terms in Q (IC.3P are already 0{^^). Thus we have to look at only at terms 
which when acting on (JC.51) do not produce extra derivatives of metric. This happens only 
when oq coefficient is left in (IC.5|) and when "free" derivatives of (IC.4|) act by pairs on the 
world function a in fIC.SP thus producing the components of metric but not its derivatives. 
One needs to know only the following structures at coinciding arguments: V a^ pff = ga^, 
VadrO' = 0, drdrff = N"^ . Up to total derivatives this gives 



(47r(T)<i/2 



JdfM 



^L^F"{q)--(k'' + 2h,k'^) + --F"{qy—- ^ '' ! '" 



Is 



Is' 



-^F"(g)(A;^ + 2fc,,F^)--^F"(g) 






(C. 



Returning back to original metric (i.e. replacing N=^Ja / s N and kij=^ys/a kij) and com- 
bining together (1C.7I) and (IC.SP one finally comes tql3 



Tr Q e-"^ 



where (i/i = drd'^'hc N^/g. 

Integration over auxiliary proper time a (which gives just certain powers of q) and subse- 
quent integration over q completes calculation of (I6.24p . (I6.25p . 



^^ We independently confirmed this result by performing a perturbative calculation using momentum repre- 
sentation (similar to that of Appendix A but involving second order in perturbation) . 
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